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Алгоритм выбора интервалов загрубления сетки моделирования, предложенный в статье, позволяет выявлять наиболее одно- 
родные по фильтрационно-емкостным характеристикам смежные слои трехмерных цифровых моделей месторождения нефти 
и газа и минимизировать тем самым размерности моделей. Авторский алгоритм наиболее эффективен для анализа и загрубле- 
ния моделей, сформированных при помощи стохастических методов, и опробован на примере моделей нефтяных месторожде- 
ний Томской области. 


Моделирование с применением современных 
информационных систем и технологий сопровожда- 
ет жизненный цикл месторождения от момента по- 
иска углеводородного сырья до ликвидации скважин 
после их отработки. Для проектирования и опера- 
тивного управления разработкой создаются трехмер- 
ные цифровые модели, на основе которых принима- 
ются управленческие решения. Такие модели имеют 
ограниченную техническими и временными ресур- 
сами размерность и должны максимально достовер- 
но отображать объект исследований и разработки. 

В процессе моделирования последовательно 
создаются геологическая (ГМ) и гидродинамиче- 
ская (ГДМ), как правило, ячеистые модели место- 
рождения. Первая из них является статической и 
характеризует фильтрационно-емкостные характе- 
ристики (ФЭС) составляющих объект разработки, 
вторая, рассчитанная на основе первой, - динами- 
ческой и описывает процессы фильтрации пласто- 
вых флюидов в процессе эксплуатации месторож- 
дения. Очевидно, что второй вид моделей требует 
значительно больших компьютерных ресурсов, по- 
скольку выполняется огромный объем вычислений. 


Для уменьшения времени расчета гидродинамиче- 
ских процессов на основе геологической модели 
используют загрубление исходной сетки моделиро- 
вания, т. е. уменьшают число ячеек сетки моделиро- 
вания путем объединения/осреднения смежных 
ячеек исходной сетки геологической модели. 

Загрубление сетки моделирования, которая, как 
правило, имеет вид регулярной сетки, именуемой 
гридом, должно вестись с учетом характера распро- 
странения в модели ФЭС, каковыми являются по- 
ристость, проницаемость, нефтенасыщенность, 
определяемые по данным сейсмических исследова- 
ний и информации из скважин. Более корректно 
объединять слои модели, поскольку геологическая 
среда носит слоистый характер и важно не нарушить 
характер гидродинамической связи между слоями. 

В применяемых на сегодняшний день для моде- 
лирования месторождений нефти и газа програм- 
мных комплексах и системах имеется ряд встроен- 
ных функций для объединения ячеек (пропорцио- 
нально по среднему арифметическому и т. и.), ко- 
торые не предлагают пользователю решения по 
выбору способа объединения. Специализирован- 
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ные модули выполняют осреднение соседних яче- 
ек, не сохраняя исходную слоистость модели, а на 
таких моделях значительно замедляются расчеты 
фильтрационных процессов [1—5]. Кроме того, при 
моделировании ФЭС по пласту применяют детер- 
минированный и стохастический подходы. В пер- 
вом случае, восстановление показателей в меж- 
скважинном пространстве происходит на основе 
набора нерегулярных данных (например, значения 
в скважинах) с применением интерполяционных 
методов. Наиболее часто применимы методы: три- 
ангуляция Делоне, сплайн-интерполяция, Кри- 
гинг, метод обратных расстояний [5]. Указанный 
подход применяется, если скважинные данные ох- 
ватывают вариации свойства модели (рис. 1, а). 

В случае, если моделируемое свойство не имеет 
ровных средних значений между скважинами, ис- 
пользуют стохастический подход и соответствую- 
щие алгоритмы, например, последовательную гаус- 
сову имитацию, которая использует вариограмму и 
распределение свойств входных данных для добав- 
ления дополнительных характеристик свойству 
(рис. 1, б). Модели свойств, созданные последова- 
тельной гауссовой имитацией, выглядят более кор- 
ректными, чем те, которые созданы с использова- 
нием методов детерминированного усреднения. 
Недостаток этого метода в том, что выбор варио- 
грамм носит субъективный характер и требует опы- 
та и развитой интуиции специалиста. Тем не менее, 
этот подход все более часто применяют при моде- 
лировании месторождений нефти и газа. 



Рис. 1 . Восстановление ФЭС пласта по данным в скважинах 
с применением: а) детерминированного и 6) стохас- 
тического подходов 


Выполнять анализ распределения ФЭС цифро- 
вых моделей, рассчитанных с применением стохас- 
тических алгоритмов, при помощи программных 
средств моделирования месторождений нефти и 
газа (компаний 8сЫшпЪег§ег, Ьапбтагк, ВЕНСИС 
и др.) для принятия решения о параметрах загру- 
бления сетки весьма затруднительно. Поэтому в 
лаборатории моделирования месторождений неф- 
ти и газа института «Кибернетический центр» Том- 
ского политехнического университета был предло- 
жен алгоритм выбора интервалов загрубления сет- 
ки моделирования при переходе от ГМ к ГДМ с 
применением двумерных парных корреляций, осу- 
ществляющий минимизацию размерности трех- 
мерных моделей нефтегазовых месторождений. 

Рассмотрим работу указанного алгоритма на 
примере трехмерной цифровой геологической мо- 
дели объекта разработки нефтяного месторожде- 
ния, построенной в среде программного обеспече- 
ния ІКАР КМ8 компании Кохаг [4]. Размерность 
исходной сетки моделирования: по ХхУ - 124x99, 
шаг 100 м; по 2 - 42 слоя. Число активных ячеек 
(участвующих в фильтрации) - 93115. Пусть требу- 
ется получить 10 слоев в затрубленной модели. 

При моделировании разработки месторожде- 
ний нефти важна выработка запасов и поэтому де- 
тальность описания должна быть наибольшей в 
нефтенасыщенных слоях, а поскольку в рассматри- 
ваемой модели число активных ячеек в нижних 
слоях 32-42 не превышает 800 и средняя нефтена- 
сыщенность в них составляет 0,2, указанные слои 
были объединены в единый слой гидродинамиче- 
ской модели, а верхние - основные нефтенасы- 
щенные слои детально исследованы на тесноту 
связи физических характеристик. 

Слои ячеистой модели - это, по сути, регуляр- 
ные сетки, характеризующие поля параметров - 
ФЭС модели. Анализируемыми параметрами явля- 
лись водонасыщенность, пористость и песчани- 
стость. Проницаемость была исключена из списка 
анализируемых параметров, поскольку построение 
указанного параметра в геологической модели ве- 
лось на основе пористости при помощи монотон- 
ной функции и, следовательно, характер измене- 
ния проницаемости в пространстве аналогична ха- 
рактеру изменения пористости. 

На рис. 2 показаны в виде регулярной сетки с 
градиентной раскраской 1 и 2 слои геологической 
модели по показателю коэффициента пористости. 
На рисунке отмечены также контуры залежи. Здесь 
и далее нумерация слоев модели ведется сверху (от 
кровли пласта) вниз (к подошве пласта). 

Для исследования степени схожести слоев гео- 
логической модели с целью последующего загру- 
бления и определения границ слоев сетки гидроди- 
намического моделирования применялся корреля- 
ционный анализ [6], а, именно, выполнялся расчет 
парного коэффициента корреляции для регуляр- 
ных сеток. 
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Рис. 2. Слои 1 ( слева) и 2 ( справа) геологической модели ( коэффициент пористости) 


На рис. 3 показан пример рассчитанной регу- 
лярной сетки, характеризующий парный коэффи- 
циент корреляции 1 и 2 слоев геологической моде- 
ли по значению пористости. 



Рис. 3. Поле коэффициента парной корреляции пористости 
слоев 1 и 2 и профиль сечения 

Аналогичные вычисления проводились для всех 
последовательно расположенных пар слоев модели 
(1 и 2 слои; 2 и 3; 3 и 4 и т. д.). 

Визуально и по статистическим показателям 
слоев модели сложно оценивать полученные зави- 
симости. Так, на рис. 3 линией показан профиль, по 
которому проводилось сечение полученных регу- 


лярных сеток парного коэффициента корреляции, 
на рис. 4 профиль по трем последовательно рассчи- 
танным регулярным сеткам (слои 1-2, 2-3, 3-4). 


Профиль по линии 



Расстояние, метры 

Рис. 4. Профиль по регулярным сеткам парного коэффици- 
ента корреляции параметров слоев 1-2, 2-3 и 3-4 

Количественными оценками степени тесноты 
можно считать показатели, именуемые условно сум- 
марными объемными показателями, и равные объе- 
му фигуры, ограниченной корреляционной поверх- 
ностью внутри контура - границы залежи. На рис. 5 
показана схема расчета указанного параметра. 

Таким образом, рассчитываются суммарные 
объемные показатели для параметров пористости, 
песчанистости и водонасыщенности, которые по- 
казаны на рис. 6 (по оси X отмечены слои модели, 


Площадь + Высота 



Ѵ=8823670,51 об. ед. 

Рис. 5. Схема количественной оценки тесноты связи параметров слоев модели 
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по оси У - разница значений объемных суммарных 
показателей). На графиках необходимо оценить ха- 
рактер изменения свойств от слоя к слою, т. е. вы- 
делить точки, в которых наблюдаются минималь- 
ные значения суммарных объемных показателей 
(т. е. границы слоев с минимальной схожестью) и, 
следовательно, будут границами слоев затрублен- 
ной модели. 

Прокоррелировав эти показатели, получаем 
график (рис. 7), характеризующий степень похоже- 
сти слоев цифровой геологической модели по 
ФЭС. Поиск 8-ми минимальных значений на гра- 
фике дает границы интервалов загрубления моде- 
ли, показанные также на рисунке. 


1500 
1000 
500 
0 

—о— пористость —о— песчанистость 

—о— водон асы ще н н ость 

Рис. 6. Г рафики значений суммарных объемных показате- 
лей парных корреляций параметров слоев модели 




Рис. 7. Г рафик корреляции объемных показателей и грани- 
цы варианта загрубления (минимизация коэффици- 
ента корреляции) 

Однако поиск минимумов по всему интервалу 
значений дает некорректный результат, т. к. воз- 
можна значительная неоднородность разбиения, 
что и отображено на рис. 8, где 5-й слой ГДМ со- 
ставил 12 слоев исходной ГМ -6 м. Поэтому был 
выполнен поиск минимумов внутри интервалов 
ограниченной длины. А именно, расчет проводил- 
ся следующим образом. 

Пусть п и ^-количество слоев ГМ и ГДМ, і - те- 
кущий слой ГДМ, к {]і и к ік - начало и конец интерва- 
ла поиска минимума, к тШ - текущий интервальный 
минимум, к оЧ — суммарный объемный показатель. 

Тогда для первого интервала: 

*оі=1; К=п/п І + к т \ ^ш п1 =/ шіп{^ обі }, гдеу-^Аі- 


Для последующих интервалов: 

к Л =К+ 1 ; к к =к п +0,1(п- к^/ (пЧ -\ ) 2 ; 
к тіп гГ тт{^}, гд &]=к и ,Ц, і=2,п'. 

Таким образом, рассчитан вариант загрубления 
модели с поинтервальной минимизацией коэффи- 
циента корреляции. На рис. 8 показаны получен- 
ные интервалы загрубления. 

Степень сходимости по запасам в затрубленной 
по различным вариантам модели (% изменения 
объема запасов нефти) приведена в таблице. На 
рис. 9 приведены полученные интервалы объеди- 
нения слоев ГМ в принятом варианте загрубления 
при построении ГДМ. 



Рис. 8. Г рафик корреляции объемных показателей и грани- 
цы варианта загрубления (поинтервальная миними- 
зация коэффициента корреляции) 




пористость 

во до н асы ще н н ость 


песчанистость 


Рис. 9. Интервалы объединения слоев ГМ в принятом вари- 
анте загрубления 


Таблица. Сравнение вариантов разбиения 


Варианты разбиений 

% от запасов по 
ГМ 

Пропорциональное 

7 

По статистическим параметрам 

5 

По минимизации коэффициента корреляции 

3 

По интервальной минимизации коэффици- 
ента корреляции 

2 


Из таблицы видно, что степень сходимости мо- 
дели по запасам повышается от варианта к вариан- 
ту и имеет разброс от 7 до 2 %. Несмотря на неболь- 
шой разброс в значениях и проценте отклонения 
по запасам от ГМ, применение предложенного 
способа загрубления геологической модели при 
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Начал 


Мг(ахЬхп) 
СЬнк, г/, і=1 



і=1,п 



М Г д/(ахЬ)=Мп(ахЬ) сог Мп+\(ахЬ) 


^/(бвнк, Мщі) 


Поиск поинтервальных 
минимумов, К т [ п (п) 



^Чпіп (}1 ) 

Конец 



МтіахЪхп ) - трехмерная цифровая 
геологическая модель 
а, Ъ, ип- число ячеек ГМ по оси 
Х,У и 2 

СЬнк - контур водонефтяного 

контакта — границы залежи 
нефти 

п - число слоев затрубленной 

гидродинамической модели 
МгдііахЪ ) и Мг!(ахЪ ) - і - й слой 
ГДМ и ГМ 

К т [ п (п) - границы интервалов 
загрубления ГМ 


Рис. 10. Алгоритм выбора интервалов загрубления сетки моделирования с применением двумерных парных корреляций 


построении ГДМ, оправдано, поскольку таким об- 
разом объединяются максимально однородные 
смежные слои, что влечет минимальную деформа- 
цию линий тока в фильтрационной модели. 

Таким образом, минимизация размерности трех- 
мерных моделей осуществляется по алгоритму выбо- 
ра интервалов загрубления сетки моделирования при 
переходе от ГМ к ГДМ с применением двумерных 
парных корреляций, который представлен на рис. 10. 

Описанный алгоритм позволяет выполнять ана- 
лиз статистических показателей фильтрационно- 
емкостных характеристик трехмерных цифровых 
месторождений нефти и газа, построенных при по- 
мощи стохастических методов. Поскольку такие 
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модели не поддаются визуальной оценке, а стати- 
стические характеристики (минимальные, макси- 
мальные, средние значения, дисперсия и т. п.), в 
том числе полученные внутри залежи не позволяют 
адекватно оценить степень однородности слоев мо- 
дели, поэтому предложенный алгоритм позволяет 
обоснованно выделять границы загрубления слоев 
трехмерных цифровых моделей месторождений. 

Описанный в данной статье алгоритм выбора 
интервалов загрубления сетки моделирования при 
переходе от геологической к гидродинамической 
модели был эффективно опробован при выполне- 
нии проектной документации по двум нефтяным 
месторождениям Томской области [7]. 
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